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ABSTRACT 

A modified gravity (MOG) theory, that has been successfully fitted to galaxy rotational velocity 
data, cluster data, the Bullet Cluster 1E0657-56 and cosmological observations, is shown to be in good 
agreement with the motion of satellite galaxies around host galaxies at distances 50-400 kpc. 
Subject headings: cosmology: theory — dark matter — galaxies: halos — galaxies: structure 



1. INTRODUCTION 

Dark matter has been used successfully to explain the 
rotational velocity data for galaxies and clusters of galax- 
ies, lensing by galaxies and clusters of galaxies, and cos- 
mological observations. Within the framework of gen- 
eral relativity (GR), this standard model of cosmology 
(ACDM) provides a successful description of the uni- 
verse. The model predicts a universe that is composed of 
~ 4% normal baryonic matter and ~ 24% non-baryonic 
cold dark m atter (CDM), while t he rest is made up of 
dark energy (jSpergel et al. I l2006h . However, dark mat- 
ter has not been detected after much experimental effort. 
It is possible that exotic dark matter does not exist in 
sufficient quantities to explain the large amount of as- 
trophysical and cosmological data. Instead, Newtonian 
gravity and GR must be modified to account for the data. 

A relativistic MOG has been proposed to expl ain the 
astro-physical and large scale cosmologica l data ([Moffat 
[20051 120061 [20071 iMoffat fc Tothl [2007al lbT) without ex- 
otic dark matter. The simplest version of the MOG 
is a scalar-tensor- vector grav ity (STVG) theory ([Moffat 
[20061 IMoffat fc Tothl l2007bf >. The theory has been fit- 
ted remarkably well to the rotatio nal velocity curve data 
of a large number of galaxie s ([Brownstein fc Moffat I 
I2006at IMoffat fc Tothl l2007bh . to the velocity disper- 
sion of globular clusters (|Moffat fc Tothl 120081) and to 
the mass and thermal profiles of clusters of galaxie s 
([Brownstein k Moffat I l2006br IMoffat fc Tothl l2007bf >. 
It h as also been fitted to the recent data o f Clowe et 
al. ([Clowe et al. 1 120061 : iBradac et al. 1I2006D describing 
the dram atic merging of two clusters 1E0657-56 (Bullet 
Cluster) ([Brownstein fc Moffatl [2007| ). A successful fit 
to the CMB acoustical power spectrum data and mass 
power spectrum data has been achieved with MOG and 
a possible explanation of the acceleration o f the universe 
appears naturally in the gravity theory ([Mo ffat 2007; 
IMoffat fc Toth~ll2007al). 

Recently. IKlvpin fc Pradal (|2007f ) presented an analy- 
sis of galaxy observations of the Sloan Digital Sky Survey 
(SDSS, http : //www. sdss . org/) to test gravity and dark 
matter in the peripheral parts of galaxies at distances 50- 
400 kpc from the centers of galaxies. This field of extra- 
galactic astronomy provides on e of the main arguments 
for the presence of dark matter (|Zaritskv fc Whitell994l : 
iPrada et al. 1 12003). In the present work, we repeat part 
of this analysis and confirm that MOG predicts satellite 



galaxy velocities that are consistent with SDSS observa- 
tions. 

2. DATA ANALYSIS 

The SDSS provides imaging data over 9500 deg 2 in 
five photometric bands. Galaxy spectra are determined 
by the CCD imaging and the SDSS 2.5m telescope 
on Apache Point, New Mexico. Over half a million 
galaxies brighter than Petrosian r-magnitude 17.77 over 
7400 deg 2 are included in the SDSS data with a redshift 
accuracy better than 30 km/s. 

We followed the analysis of IKlvpin fc Prada I (|2007l ) 
and obtained a data set of nearly 700000 galaxies 
from the latest SDSS da t a rele ase (Data Release 6, 
lAdelman-McCarthv et al. I (|2008[ )). We checked our 
query results by searching for, and locating several galax- 
ies the parameters of which are known from the litera- 
ture. 

Rest frame absolute magnitudes can be computed 
in the SDSS g-band from the extinction-corrected ap- 
parent magnitudes assuming a Hubble constant Hq ~ 
71 kms^Mpc -1 . All magnitudes were fc-corrected to 
z = 0. SDSS redshifts are heliocentric; before calculat- 
ing distances, we also took into account the solar system's 
motion relative to the CMB. 

Due to the complications of calculating mod i fied g ray- 
ity for non-spherical obiects. lKlvpin fc Pradal (|2007l ) re- 
stricted their analysis only to red galaxies, the vast ma- 
jority of which are cither elliptical galaxies or are domi- 
nated by bulges. We followed a similar strategy, restrict- 
ing our selection of candidat e host galaxies to luminous 
red ga laxies, as identified bv lAdelman-McCarthv et al. I 
( 2008) . Wc also restricted our selection to galaxies with 
a recession velocity between 3000 km/s and 25000 km/s, 
which yielded approximately 234000 galaxies in total. 

Still continuing wit h the strategy laid out by 
IKlvpin fc Pradal (l200l . we identified a host galaxy as 
a galaxy that is at least 4 times brighter than any other 
galaxy within a projected distance of R = 1 Mpc and 
line-of-sight velocity within 1500 km/s that of the candi- 
date host. Thereafter, we identified any galaxies at least 
4 times fainter than the candidate host within a pro- 
jected distance of 1 Mpc and with a line-of-sight velocity 
within 1500 km/s of that of the host as a candidate satel- 
lite galaxy. This strategy yielded a total of ~3600 host 
galaxies with ^8200 satellites. The procedure assigned a 
small percentage (~ 1.5%) of satellites to multiple hosts; 
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Fig. 1. — Line-of-sight velocities (top left) and number densities (bottom left) of candidate satellite galaxies as a function of projected 
distance from the candidate host. After removal of candidate interlopers, the line-of-sight velocities (top right) and number densities 
(bottom right) are also shown. The solid lines in the top plots are moving window averages. 



we made no attempt to eliminate these from the data 
set. 

The velocity distribution of candidate satellites as a 
function of projected distance from candidate hosts is 
shown in Figure[T] (top left). We note that our choices of 
parameters (in particular, our choices of the factor four 
in the difference in luminosity between candidate hosts 
and satellites, or the restriction of hosts to luminous red 
galaxies) did not substantially affect the resulting distri- 
bution. 

Obviously, the procedure outlined here does not guar- 
antee that all candidate satellites are, in fact, true satel- 
lites of their respective candidate hosts. Indeed, one 
would expect a mix of satellite galaxies and a random 
background of interlopers. The number of interlopers 
should increase as a linear function of the projected dis- 
tance, whereas the number of true satellites should peak 
at some typical satellite orbital radius. This expectation 
is confirmed if we plot the number density of candidate 
satellites as a function of projected distance (Figure Q] 
bottom left). We are immediately presented with a pos- 
sible editing strategy: after binning the data according 
to projected distance, we can remove a number of can- 
didate satellites from the sample to account for the ran- 
dom background. By assuming that at 1 Mpc, none of 
the candidate satellites are, in fact, true satellites (i.e., 
all are random interlopers), we obtain a simple formu- 
lation: the number of galaxies to be removed from each 



bin is -ZViMpc-R, where -/Vim P c is the number of candidate 
satellites at 1 Mpc. The resulting number density of can- 
didate satellites is shown in Figure [1] (bottom right). 

One question remains: what is a satellite? In other 
words, which galaxies are to be removed from the sam- 
ple? We answer this question by noting that at a large 
projected distance from a candidate host, a candidate 
satellite is much less likely to be a true satellite of that 
host if its line-of-sight velocity also differs greatly from 
that of the host. Therefore, the candidate satellites that 
we remove from the sample are those with the greatest 
difference in line-of-sight velocity relative to the respec- 
tive candidate host. The velocity distribution of the re- 
maining satellites, shown in Figure [1] (top right), is the 
data set that we study. As this simplistic selection pro- 
cedure tends to underestimate satellite velocities at high 
projected distance, in the following, we shall restrict our 
investigation to a projected distance of 400 kpc or less. 

3. MODIFIED GRAVITY PREDICTIONS 

Predictions for modified gravity can be made by solv- 
ing the Jeans equation, which gives the radial veloc- 
ity dispersion of ( r) as a function of radial distance. 
iKlvpin k. Prada 1 (|2007| ) find that neither Newtonian 
gravity, nor MOND are compatible with observations. 
lAngus et al. I ()2008| ). however, demonstrated that a suit- 
ably chosen anisotropic model and appropriately cho- 
sen galaxy masses can be used to achieve a good fit for 
MOND. 
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Radial velocity dispersions in a spherically symmet- 
ric gravit ational field can be comput ed using the Jeans 
equation (B innev fe Tremaine. 1119871 ): 



d(va A r 
dr 



2v 



d<Z> 



(1) 



where v is the spatial number density of particles, v r is 
the radial velocity, (3(r) — 1 — [a 2 (r)+a 2 > (r)]/2a 2 (r) is the 
velocity anisotropy, $(r) is the gravitational potential, 
and we are using spherical coordinates r, ef>. We can 
write Eq.|T|) in the form 



Here, we have 



da 2 Aa 2 r 
dr r 



A = 2/3(r)+~/(r), 



(2) 
(3) 



where "/(r) = d\nv(r) / dhir. 

If we assume that the velocity distribution of satellite 
galaxies is isotropic, (3 = 0. In general, j3 needs to be nei- 
ther zero nor constant. The number density of candidate 
satellites favors a value of 7 ~ —2.5. We note that our 
MOG predictions are not particularly sensitive to these 
parameters; no "fine-tuning" is needed to obtain good 
agreement with the data. 

The observed velocity dispersion is along the observer's 
line of sight, seen as a function of the projected distance 
from the host galaxy. Therefore, it is necessary to inte- 
grate velocities along the line-of-sight: 



/ [y 2 + (1 - (3)R 2 ] r- 2 o 2 Mv{y) dy 



(4) 



/ v{v) dy 




where v is the spatial number density of satellite galaxies 
as a function of distance from the host galaxy, and y is 
related to the projected distance R and 3-dimensional 
distance r by 

V 2 - (5) 



r 2 = R 2 



Changing integration variables to eliminate y, we can 
express the observed line-of-sight velocity dispersion as 
a function of projected distance as 



J{r 2 - f3R 2 )a 2 {r)v{r)/r^r 2 - R 2 dr 
3os(*) = - 58 ■ (6) 



J rv(r)/y/ r 2 
R 



R 2 dr 



From the field equations derived from the STVG 
action, we obtain the modified Newto nian ac c elera- 
tion law for weak g ravitational fields (|Moffat 2006; 
iMoffat fc Toth Il2007bh of a point source with mass M: 



SmogM 



G N M 



{l + a[l-e-^(l+/xr)]}, (7) 



where Gat is the Newtonian gravitational constant, while 
the MOG parameters a and fi determine the coupling 
strength of the "fifth force" vector M to baryon matter 
and the range of t he force, respectively. 

In recent work (jMoffat fc Toth Il2007bt) . we have been 
able to develop formulae that predict the values of the a 
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Fig. 2. — Root-mean-square line-of-sight velocities and theoreti- 
cal predictions. Line of sight velocities are binned by projected dis- 
tance, in bins of 40 kpc. The parameter-free MOG prediction (solid 
blue line) using nominal host galaxy masses is in good agreement 
with the data. MOND (dotted red line) is in marginal agreement, 
which can be improved by, for instance, varying the parameters (3 
and 7 with radial distance. The same is not true for the Newto- 
nian prediction (dashed green line). The two panels show two host 



galaxy luminosity ranges. 
-21.1 > Ml > -21.6. 
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and /j, parameters, in the form 
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(8) 
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(12) 



are determined from galaxy rotation curves and cosmo- 
logical observations (|Moffat fc Toth H2007bT) . 

The MOND acceleration (?mond is given by the solu- 
tion of the non-linear equation 



.9MONDM 



|#mond|\ _ G N M(r) 



no 



where M is the mass of only baryons and a = 
The form of the function _u(x) 



(13) 

1.2 x 
origi- 
nally proposed by Milgrom (|Milgrom I fl983h is given 
by /j,(x) = x/y/l + x 2 ; however, better fits and better 
asy mptotic behavior are a chieved using )J,(x) = x/(l + 
x) (|Klypin fc Pradall2007l ). 

Following in the footsteps of lKlvpin fc Pradal (|2007f ). 
we grouped satellite galaxy velocities for host galaxies 
in two luminosity ranges: —20.5 > Mg > —21.1, and 
—21.1 > Mg > —21.6. The corresponding masses for the 
host galaxies, calculate d by iKlypin fc Prada I (|2007l ) on 
the basis of the work of lBell fc de JongT ( 200lD . are 



Ml = 7.2 x 10 10 M (• 
1.5 x 10 11 M Q (- 



M 2 



-20.5 > Ml > 
-21.1 > M 2 > 



-21.1), 
-21.6). 



(14) 
(15) 
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We plotted the MOG prediction using these nominal 
masses and the values of (3 — 0, 7 = —2.5, along with 
the predictions of MOND and Newtonian gravity using 
the same parameters. As Figure [2] shows, the MOG pre- 
diction is consistent with observational data, although 
MOG prefers a somewhat higher than nominal mass for 
the hosts; a least-squares fit yields M* ~ 9.3 x 10 10 M 
and Ml ~ 2.2 x 10 11 M . 

4. CONCLUSIONS 

Observational data presented by 

Adelm an-McC arthv et al. I (|2008t) and studied by 



Klvpin fc Pradal (|2007fl are viewed as evidence of the 



success of the ACDM model. However, the data are 
also in excellent agreement with the predictions of the 
parameter-free solutions of our modified gravity theory, 
without postulating a peculiarly distributed halo of 
exotic dark matter. 
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